Measuring the geometry of the Universe in the presence of isocurvature modes 
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The Cosmic Microwave Background (CMB) anisotropy constrains the geometry of the Universe 
because the positions of the acoustic peaks of the angular power spectrum depend strongly on the 
curvature of underlying three-dimensional space. In this Letter we exploit current observations to 
determine the spatial geometry of the Universe in the presence of isocurvature modes. Previous 
analyses have always assumed that the cosmological perturbations were initially adiabatic. A priori 
one might expect that allowing additional isocurvature modes would substantially degrade the con- 
straints on the curvature of the Universe. We find, however, that if one considers additional data 
sets, the geometry remains well constrained. When the most general isocurvature perturbation is 
allowed, the CMB alone can only poorly constrain the geometry to Qo ~ 1-6 ± 0.3. Including large- 
scale structure (LSS) data one obtains fio = 1-07 ± 0.03, and Qq — 1.06 ± 0.02 when supplemented 
by the Hubble Space Telescope (HST) Key Project determination of Hq and SNIa data. 



One of the most striking successes of observational cos- 
mology over the past few years has been the constraint 
of possible departures of the geometry of our universe 
from a spatially flat (Euclidean) geometry. Such a de- 
parture from spatial flatness may be characterized by the 
dimensionless quantity ^k, related by the Friedmann- 
Robertson- Walker equations to the fractional energy den- 
sity (with respect to the critical density) fig of the Uni- 
verse by the relation fix = 1 — Oq. If < 0, the 
Universe has a three-dimensional spherical geometry, if 
> 0, a hyperbohc geometry, and if = 0, a Eu- 
clidean (flat) geometry. When a structure of a given 
physical size is viewed from a cosmological distance, its 
angular size depends sensitively on the spatial geometry. 
Different spatial geometries yield different laws of per- 
spective. Such differences may be exploited to obtain a 
measurement of Q,k using the angular power spectrum of 
the CMB anisotropy Q . In the standard Big Bang ther- 
mal history, the Universe, previously ionized, recombined 
into neutral atoms at a redshift Zrec ~ 1080. At this mo- 
ment the photon mean free path increased precipitously, 
and since then the CMB photons have travelled freely 
toward us with virtually no rescattering. Prior to re- 
combination the previously tightly coupled plasma had 
undergone acoustic oscillations, manifested in the CMB 
spectrum today by the so-called Doppler peaks. At re- 
combination these oscillations were characterized by a 
physical scale d*^^, the size of the sound horizon at re- 
combination. In order to convert d*^^ into an angular 
scale 9*, another length scale dAD is required. Let dLs 
be the diameter of the last scattering surface expressed 
in terms of present day co-moving units. In a flat geom- 
etry, it follows that 9* = {zrec + l)c?*ec/c'LS- Howcvcr, 
in a curved geometry, dLS must be replaced with the 
apparent angular diameter distance d^D, defined as 
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to account for the (de)focusing of rays by the nonzero 
spatial curvature 

The detection of the acoustic peak in the CMB by 
a number of experiments Q and more recently with the 
Wilkinson Microwave Anisotropy Probe (WMAP) "T] has 
been used to constrain the spatial curvature of the Uni- 
verse. These analyses, however, assumed that the pri- 
mordial perturbations were adiabatic — that is, at very 
early times the universe was governed by a common, spa- 
tially uniform equation of state and all the components 
contributing to the stress-energy shared a common pecu- 
liar velocity field. When this assumption is relaxed other 
modes, the so-called isocurvature modes, arise, and their 
presence generically alters the positions of the Doppler 
peaks used to determine the spatial curvature. It is of 
interest to learn whether the constraints on Qk weaken 
significantly when the assumption of adiabaticity is re- 
laxed. Previous studies have projected that without as- 
suming adiabaticity the ability to determine ^Ik with 
the WMAP data would be severely compromised In 
this Letter we examine this question and find that 
is poorly constrained by the CMB data alone; however, 
when LSS structure data is included, one is able to im- 
pose stringent constraints on ^k- 

The data sets considered here include the CMB 
anisotropy data and the SDSS determination of the 
galaxy power spectrum, as well as the determination of 
the Hubble constant and the luminosity-redshift relation 
for Type la Supernovae. The fluctuations of the CMB are 
characterized by the angular power spectrum C'i where 
(|AT/„p) = £{£ + l)Ce/2TT, which may be transformed 
into an angle using i « tt/9. For a cosmology with adia- 
batic initial cosmological perturbations (and other initial 
conditions as well), the angular CMB power spectrum C; 
exhibits a series of peaks and troughs whose positions 
scale according to the above determination of 9*. For 
a Euclidean universe of an age consistent with current 
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FIG. 1: The CMB temperature spectra for the adiabatic and 
isocurvature modes AD, Cl, nid, niv for fio = 0.8,1.0,1.2, 
where cot = 0.022, coa = 0.12, r = 0.1, = 1, /i = 0.7. 

measurements of the Hubble constant and adiabatic 
initial conditions, the first so-called Doppler peak is sit- 
uated at e fa 220, as observed by WMAP [3|, leading 
to constraints, assuming adiabatic initial conditions of 
no = 1.09±!; ?§ 0, or 0.98 < Qq < 1.08 at 95% confi- 



The 



dence, assuming a Hubble constant of h > 0.5 
eometry was further constrained to ^Iq = 1.02 ± 0.02 
~\ by including data from the 2dFGRS or SDSS galaxy 
surveys |^, Il0j |. supernovae redshift- luminosity relation 
[Tl| and the HST key project determination of /i 
We now consider more general initial conditions where 
isocurvature modes and possible correlations among 
themselves and with the adiabatic mode are also al- 
lowed. Under this framework four additional perturba- 
tion modes are possible: the CDM (CI), baryon (BI), 
neutrino density (NID) and neutrino velocity (NIV) 
isocurvature modes, discussed in detail in 0, llJ. Con- 
straints on these modes have been presented in |l3l 
for flat cosmologies. Each of these modes predicts distin- 
guishable CMB power spectra templates shown in Fig. 
n except for the BI and CI modes which have almost 
identical spectra. An admixture of new modes can shift 
the peak positions and weaken the curvature constraint. 
We consider a family of cosmological models character- 
ized by five free parameters, the curvature f2fe, the baryon 
density parameter uih — flhh^, the cold dark matter den- 
sity uJd = n.dh?' , a cosmological constant f^A and the op- 
tical depth to recombination t subject to the constraint 
r < 0.3, all which are assigned uniform priors. We impose 
the weak priors < fib < 1 and < fi^ < 3 and con- 
sider neither a varying dark energy equation of state nor 
neutrino masses. With four distinguishable scalar per- 
turbation modes, labelled by i, j = AD, CI, NID, NIV, the 
initial perturbation power spectrum consists of ten sym- 
metric matrix elements A*-' {k), parametrized as A*-' (k) oc 

ik/kf)) 2 1^ introducing a scalar spectral index Ui for 
each mode and a pivot scale taken as fco — 0.05/i/Mpc. 
The angular power spectrum of temperature and polar- 
ization anisotropics corresponds to the quadratic observ- 
able 



rik 

-A^^{k)eiAk)eiBik) (2) 




0.0 0.2 0.4 0.6 0.0 0.2 0.4 0.6 

FIG. 2: Marginalized 68% (pale grey) and 95% (dark grey) 
constraints on Qrn and Ha for adiabatic (left, AD) and mixed 
initial condition models (right, AD+ISO). The top two rows 
use CMB data, the bottom row includes LSS data {h > 0.5) 
and shows constraints from SN la data alone (dotted). 

where 0£^(fc) is the radiation transfer function for the 
mode i and observable A — T or E, T indicating the 
temperature anisotropy and E the electric polarization. 
We compute the ten components of the theoretical CMB 
temperature and polarization spectra and matter 
power spectra (fc) for each correlation using the pub- 
licly available CAMB package ^3 • To generate all ten sets 
of spectra for / < 1600 and Pm{k) up to fc = 0.7/i/MPc 
takes about 20 sec for a flat model and 70-100 sec for 
curved models using two 3GHz processors. We modified 
the code to sum the spectra with varying amplitudes us- 
ing the methods described in 14], including nine relative 
amplitudes and one overall amplitude (or one amplitude 
for adiabatic models). The ten coefficients quantify 
the relative contribution to the power in the CMB by 
each correlation, with X^i j=i '^fj — 1- -^ad has j, j — AD. 

We compare models to the WMAP temperature and 
temperature-polarization CMB data covering 2 < £ < 
800 4] using the likelihood function in 16], together with 
small-scale CMB data covering 800 <1 < 1600 [Hill. 
To increase computational speed we do not use £ > 1600 
data. We also add the SDSS galaxy power spectrum 
data Pg{k) in the linear regime k < 0.2VMpc llO|. We 
include a bias parameter b such that Pg{k) — b^Pm jk) 
and use the non- linear approximation for Pm{k) in [I3 
implemented in CAMB. Finally, the analysis is repeated 
including the likelihood as a function of 0^ and f^A in- 
ferred from SN la redshift-luminosity measurements [ll| , 
and also the Gaussian prior on h of 0.72 ± 0.08 obtained 
by HST JJ. The posterior distributions are sampled us- 
ing the Markov Chain Monte Carlo methods described in 
l20l . We derive fig = 1 — and the isocurvature frac- 



tion fiso = Ziso/{ziso + Zad), where = 1- 
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marginalized median and 68% confidence intervals. We 
test the software by reproducing results for flat models 
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FIG. 3: Constraints on the Hubble constant h and the frac- 
tional energy density f2o, with panels as in Fig. |5| 



Dataset 



h 







> 0.2 1.2±!j;? 

>0.5 1.03 ±0.03 

>0.2 1.04 ±0.03 

>0.5 1.04 ±0.02 

0.02 



CMB 
CMB 

CMB ± LSS 
CMB + LSS 

CMB + LSS + SNia 0.72 ± 0.08 i.oi±i;:oi 

TABLE I: Median and 68% confidence intervals for the frac- 
tional energy density flo for adiabatic models with minimum 
Hubble constant h or Gaussian HST prior 0.72 ± 0.08. 

obtained with the DASh code [Tll2ll|. 

We first explored the determination of the spatial cur- 
vature allowing only the adiabatic models, characterised 
by seven parameters (eight with LSS), and including var- 
ious combinations of data sets. With only the CMB 
data and the weak prior 0.2 < ft, < 1, the curvature 
is poorly constrained as shown in Table The CMB 
provides a good measure of the angular diameter dis- 
tance dAD — {zrec + l)c^*ec/^*i t)ut a dcgcncracy exists 
between the curvature and the conformal distance to last- 
scattering dLs 22], as illustrated in the top-left panels of 
Figs, m and |3| d^s increases with decreasing h or JIa 
and with increasing flm- The constraint on the geome- 
try is considerably tightened when h > 0.5 is imposed. 
The degeneracy is also broken by the LSS data. The 
shape of the power spectrum constrains flmh, ruling out 
the closed models with high ri„j allowed by the CMB 
data, which instead constrains Q,mh? ■ Similarly super- 
nova luminosity-redshift measurements exclude models 
with high and low JIa, as shown in Fig. |2J By in- 
cluding CMB, LSS, SNIa and HST data combined, the 
geometry is strongly constrained to ilo = 1.01±g;oi, in 
agreement with [3, Ig] ■ We now consider how the mea- 
surements are affected when the full range of correlated 
isocurvature and adiabatic initial conditions are admit- 
ted. Initially the spectral indices of the modes are fixed 
as a single parameter, giving rise to a family of models 
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FIG. 4: CMB temperature (top) and galaxy power spec- 
tra (bottom) for a curved cosmological model dominated by 
isocurvature with fiso = 90%, f2o = 1-9. It fits current CMB 
data well (TT and TE), but does not fit LSS data and is incon- 
sistent with BBN and HST observations {ujb = 0.08, h = 0.3). 
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TABLE II: Median and 68% confidence intervals for fio and 
median isocurvature fraction fiso for models with mixed ini- 
tial conditions. 

described by sixteen parameters (seventeen with LSS). 
Figs, m 13 and [3 and Table ITU demonstrate that the deter- 
mination of the geometry of the Universe is significantly 
degraded by the inclusion of these additional degrees of 
freedom. As in the adiabatic cases, where flo = 1.2±Q f 
results from the CMB data alone, the data favor a closed 
Universe with JIq = 1-6 ± 0.3 when isocurvature is al- 
lowed, with the priors h > 0.2 and Slt^ < 3 constraints 
cutting off larger curvatures. The extremely closed mod- 
els, such as the one in Fig^ are dominated by up to 90% 
isocurvature (at 2cr), and exhibit low h, high Qm and low 
flA, extending the degenerate direction observed for the 
adiabatic models. Many models dominated by the NIV 
mode require a high baryon content (wf, « 0.06), so the 
acoustic peak positions are shifted to smaller scales by 
lowering the sound speed at recombination. This effect 
contributes to the exclusion of isocurvature-dominated 
extreme open models. Closed models dominated by the 
NID or CI modes do not fit the CMB data, particularly 
at large scales. 

Additional data constraining d^s improve the curva- 
ture constraints. The prior ft > 0.5 alone suffices to ex- 
clude highly closed models. Alternatively, the LSS data 
constrains the allowed range to ^Iq = 1.07 ± 0.04 in two 
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FIG. 5: Marginalized posterior distributions for Qo for adia- 
batic models (AD,unshaded) and mixed initial conditions with 
common spectral indices (ad+iso, shaded). The dark shaded 
and dashed curves show results for h > 0.5, otherwise h > 0.2. 

ways. Firstly, its shape constrains and rules out 

high ri,„ models, since in most models the matter power 
spectrum is dominated by the adiabatic mode. Secondly, 
the LSS data rules out models dominated by NIV, which 
predict significant baryon oscillations. The supernova 
and HST data further constrain d^s- The combination 
of all datasets gives Hq = 1.06 ± 0.02, with {fiso) = 0.54. 
The distribution for fiso, however, is very broad and does 
not appear to disfavor the adiabatic models. Within this 
class of models the data prefer higher baryonic densities 
{ivt = 0.04 ± 0.02 for CMB+LSS data) than consistent 
with nucleosynthesis measurements. We find that adding 
the Gaussian prior uji, = 0.022 ± 0.002 from recent BBN 
measurements to the CMB+LSS analysis with the 
prior h > 0.5 reduces the isocurvature contribution by 
just under a half and yields ilo = 1.06 ±0.03. We also re- 
peated the analysis for an enlarged space of models where 
each of the isocurvature modes is given an independent 
tilt within the range < rii < 2. The resulting twenty- 



dimensional model constrained by CMB-)-LSS with the 
h > 0.5 prior gives fio = 1-08 ± 0.03 and the follow- 
ing constraints on the spectral indices: h-^d — 

1.1 ±0.2, 

nci — 1-1 i 0.5, tinid = 1-4 ± 0.3, n^iv — 1-1 ± 0.3, with 
isocurvature fraction (fiso) = 0.5. The region of high 
probability density, however, includes models with vir- 
tually no isocurvature. We do not find that the models 
with isocurvature modes offer a better fit than the adia- 
batic models, beyond what may generically be expected 
as a result of enlarging the number of degrees of freedom 
of the theoretical model. This result is analogous to the 
flat case considered in jl3 |. 

In this letter, we have substantially strengthened the ev- 
idence for an almost flat Universe, with the data suggest- 
ing a slight preference for a mildly closed Universe. In 
considering a very general class of initial conditions, we 
have removed the assumption of adiabaticity and hence 
on specific models of the early universe such as inflation. 
In principle this could have completely removed the abil- 
ity to constrain the geometryof the Universe with current 
data sets, as was shown in [3. Indeed we find with this 
wider range of parameters a behavior similar to what has 
already been noted for purely adiabatic initial conditions: 
that the CMB alone does not suffice to constrain the ge- 
ometry of the Universe. Once one includes additional 
information (constraining h or Qmh), it is possible to pin 
down the geometry with appreciable precision. Conse- 
quently, for very general initial conditions where pertur- 
bations are imprinted at early times and subsequently 
evolved through gravitational collapse, we conclude that 
the Universe is very nearly fiat. 
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